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ABSTRACT 

We propose a simple and unambiguous way to deduce the parameters of black holes 
which may reside in AGNs and some types of X-ray binaries. The black-hole mass and 
angular momentum are determined in physical units. The method is applicable to the 
sources with periodic components of variability, provided one can assume the following: 
(i) Variability is due to a star or a stellar-mass compact object orbiting the central black 
hole and passing periodically through an equatorial accretion disk (variability time- 
scale is given by the orbital period), (ii) The star orbits almost freely, deviation of its 
trajectory due to passages through the disk being very weak (secular); the effect of the 
star on the disk, on the other hand, is strong enough to yield observable photometric 
and spectroscopic features, (iii) The gravitational field within the nucleus is that of the 
(Kerr) black hole, the star and the disk contribute negligibly. 

Subject headings: black-hole physics — galaxies: nuclei — X-ray binaries — line: profiles 
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1. Introduction 

There is an ever-growing theoretical and observa- 
tional evidence that the dark masses present in some 
galactic nuclei and X-ray binaries are black holes 
(Rees 1998). Several ways have been suggested of how 
to deduce the properties of these black holes. The first 
rough estimates of their mass were based on energy 
considerations and limits implied by the shortest vari- 
ability timescale (e.g. Begelman, Blandford, & Rees 
1984). More precise methods became possible with 
modern observational techniques like HST, VLBI and 
X-ray satellites. The aim of the present paper is to 
propose and further discuss a method to determine 
the parameters of a black hole in a system with two 
periodic components in the observed signal which are 
due to orbital motion of a stellar object and its pas- 
sages through the accretion disk of the hole. We start 
with a brief summary of the subject. 

In galactic nuclei, the existence and size of the dark 
central mass are deduced from the curve of the cen- 
tral surface brightness of the nuclei and from the spa- 
tial distribution and dynamics of surrounding gas and 
stars (see Kormendy & Richstone 1995 for a review 
and references). In particular, the standard model 
for active galactic nuclei (AGNs) with a black hole 
and an accretion disk (e.g. Rees 1984; Blandford & 
Rees 1992) has been further supported when the nu- 
clear gas in several active galaxies was found to be in 
a gravity-dominated nearly circular motion in a disk 
(e.g. Jaffe ct al. 1993; Gallimore, Baum, & O'Dea 
1997). Rotation curves of the observed nuclear gas 
disks indicate a central compact dark mass of some 
(2.0-3.5) xlO^M© in M87 (NGC 4486) (Merritt & Oh 
1997; Macchetto et al. 1997), of about 4.9xlO^M0 in 
NGC 4261 (Ferrarese, Ford, & Jaffe 1996), and of 
3.6xlO^M0 in NGC 4258 (Miyoshi et al. 1995; Maoz 
1995). 

The observed excess of quasars at high redshifts 
has stimulated a search for black holes in quiescent 
galaxies, since they may have been both the engine 
and then the residue of a former activity (Haehnelt & 
Rees 1993). In this case the most reliable present data 
come from the observations of stellar motion (van den 
Bosch & de Zeeuw 1996). Velocity profiles of several 
stellar nuclear disks offer strong evidence for black 
holes — of 2.6xlO^M© in our own Galaxy (Genzel 
et al. 1997; cf. Mezger, Duschl, & Zylka 1996), of 
2xlO^M0 in NGC 3115 (Kormendy et al. 1996a), of 
IxlO^Mo in NGC 4594 (the Sombrero galaxy) (Ko- 



rmendy et al. 1996b), of 3xlO^'ilf3 in M32 (van den 
Marel et al. 1997), and probably also of 6x IO^Mq in 
NGC 4486B (Kormendy et al. 1997). 

However, current optical and radio studies still 
probe only regions above 10^^-=- 10^ gravitational radii 
of the putative holes which is by far not enough to 
resolve any imprints of the relativistic effects due to 
the collapsed centre, and in particular to deduce the 
rate of its rotation. 

Information from within the region of a few tens of 
gravitational radii have been brought by X-ray satel- 
lite observations of a sample of Seyfert 1 galaxies 
(Nandra et al. 1997 and references therein; see also 
Fabian 1997 for a survey, and Tanaka et al. 1995 for 
the best-quality measurement from MCG-6-30-15). 
These observations yielded the profiles of the fluo- 
rescent iron Ka emission line which is likely to be in- 
duced by illumination of the very inner parts of an ac- 
cretion disk by a source off the disk plane (e.g. Matt, 
Perola, Piro, & Stella 1992; Petrucci & Henri 1997). 
Indeed, the observed broad profile skewed to lower 
energies is best explained in terms of the Doppler 
and gravitational redshifts (Fabian et al. 1995); the 
exact origin of the illuminating source remains un- 
clear. Evidence for smallness of the emitting regions 
also comes from the resolved rapid variability of the 
Ka lines (Nandra et al. 1997, and references therein). 
Note that recently Bao et al. (1997) reminded another 
unique (and measurable) feature of black hole sources 
involving accretion disk: the energy-dependent vari- 
ability of polarization of their X-rays, originally dis- 
cussed by Stark & Connors (1977). 

A way to estimate parameters of the central black 
hole is strongly related to the source variability. It 
is assumed, in fact, that a time delay between vari- 
ations of the emission line strength of the AGN and 
that of the continuum can be due to the time travel 
between the central source and the surrounding line- 
emitting gas. Radial distance and velocity of propaga- 
tion of the disturbances impose limits on the black- 
hole mass (Blandford & McKee 1982; Krolik et al. 
1991), but the resulting estimates inherit uncertain- 
ties of the underlying assumptions. Alternatively, it 
has been proposed to deduce the masses of the central 
nuclear bodies from temporal changes of the observed 
emission lines (Stella 1990), while rotational param- 
eters from the position, intensity, width (Hameury, 
Marck, & Pelat 1994; Martocchia & Matt 1996), and 
profile of the lines (Bromley, Chen, & Miller 1997; 
Dabrowski et al. 1997; Reynolds & Begelman 1997; 
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Rybicky & Bromley 1998; cf. also Bao, Hadrava, & 
0stgaard 1994 and 1996). The X-ray photometric 
light-curve profile from a hot spot orbiting not far 
from the horizon also reflects rotation of the black 
hole (Asaoka 1989; Karas 1996). 

Similar research has been focused on galactic X- 
ray binaries where a central engine was also proposed, 
namely an accretion disk fed by overflow from a sec- 
ondary star. Here the evidence for a black hole in- 
volves non-stellar appearance ("invisibility") of one 
of the components, the presence of X-radiation, vari- 
ability, spectral features (mainly the presence of rela- 
tively strong ultrasoft component and of hard X-ray 
tail) , and of course a large lower limit on mass of the 
dark component (> 3Mq), deduced from the orbital 
parameters of the (stellar) companion (namely from 
the mass function of the system). (See the surveys 
in Lewin, van Paradijs, & van den Heuvel 1995; the 
reviews of black-hole binaries with a thorough list of 
references have also been given by Barret, McClin- 
tock, & Grindlay 1996; Charles 1997. Later results 
are due to Beekman et al. 1996 and 1997; Orosz & 
Bailyn 1997.) 

It has been proposed recently to infer the parame- 
ters of the dark component from observable behaviour 
of the disk. An observable modulation of the X-ray 
emission may be caused by disk vibrations whose low- 
est frequency depends on the mass and rotation of 
the hole (Nowak et al. 1997 and references therein). 
Zhang et al. (1997) proposed, also on the basis of 
the standard thin accretion disk model, how the an- 
gular momentum (spin) of the black hole in an X- 
ray binary could be inferred from the strength of an 
ultrasoft X-component. Very recently Cui, Zhang, 
& Chen (1998) have proposed "that certain types 
of quasi-periodic oscillations (QPO) observed in the 
light-curves of black-hole binaries are produced by X- 
ray modulation at the precession frequency of accre- 
tion disks, because of relativistic dragging of inertial 
frames around spinning black holes" . Given the mass 
of the hole, they arc able to derive its spin by com- 
paring the computed disk precession frequency with 
that of the observed quasi-periodic oscillations. This 
mechanism requires a warped accretion disk, the as- 
sumption that clearly calls for further investigation 
(Markovic & Lamb 1998). Scheme for the source 
variability discussed in the present paper is different 
from that which has been proposed for QPOs in the 
above-mentioned papers. Here, an orbiting compan- 
ion of the black hole is involved — the assumption 



which restricts variability time-scales relevant for the 
model. 

Narayan, McClintock, & Yi (1996), and Narayan, 
Garcia, & McClintock (1997) argued that the black- 
hole X-ray binaries could also be recognized accord- 
ing to a larger variation in luminosity between their 
bright and faint states than is expected in the sources 
with neutron stars. This way of identifying black 
holes follows from the low radiative efficiency of ad- 
vection-dominated accretion flows which could oc- 
cur around black holes (with greater "sucking-in" 
power and no solid surface) rather than around other 
compact objects. Various under-luminous accretion- 
powered astrophysical systems, mainly quiescent tran- 
sient X-ray sources and low-luminosity galactic nuclei, 
could be interpreted as black holes with advection- 
dominated accretion disks (see Lasota & Abramowicz 
1997 for a survey); Reynolds et al. (1996) suggested 
that advection-dominated mode of the final stages of 
accretion could also account for the "quiescence" of 
the black hole in M 87. The lack of a hard surface of 
a black hole also plays crucial role in an alternative 
explanation by King, Kolb, & Szuszkiewicz (1997) of 
X-ray transients which are dynamical black-hole can- 
didates, namely in terms of a weaker (stabilizing) ir- 
radiation of the disk by the central accreting source. 

In the present paper we propose a method appli- 
cable, under certain assumptions, to the black- hole 
sources with periodic components of variability. It 
provides, unambiguously, both the mass and the spe- 
cific angular momentum of the black hole in physical 
units. We start from the model of possible periodic 
variability of black- hole sources which considers a thin 
accretion disk lying in the hole's equatorial plane, and 
an object (a "star" ) which intersects the disk period- 
ically while orbiting about the centre.^ Such a sys- 
tem is characterized by two angular frequencies — 
that of azimuthal revolution of the star and that of 
its latitudinal oscillation about the equatorial plane. 
Both frequencies are in principle measurable at infin- 
ity, the azinuithal one from spectrophotometry, while 
the latitudinal one from photometry provided that the 
passages of the star produce a strong enough modula- 
tion of the source. Karas & Vokrouhlicky (1994) illus- 
trated, by Fourier analysis of simulated photometric 



^See e.g. Krolik et al. 1991; Mineshige, Ouchi & Nishi- 
mori 1994; Ipser 1994; Zakharov (1994); Bao, Hadrava, & 
0stgaard (1994); Kanetake, Takeuti, & Fukue (1995), and ref- 
erences therein for alternative explanations of periodic or quasi- 
periodic variability. 
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data, how tlic' respective two peaks can be recognized 
in the power spectrum. This model was cultivated no- 
tably in connection with the NGC 6814 galaxy whose 
putative variability, however, was later recognized as 
being due to a source in our Galaxy. Other possible 
targets are discussed in the current literature. For 
example, optical outbursts in the blazar OJ 287 have 
recently been modelled in terms of a black-hole bi- 
nary system by Villata et al. (1998). (The source 
exhibits several time-scales: fcaturc-lcss short-term 
variability, 12-yr cycle, and, possibly, a 60-yr cycle.) 
These authors, however, presume both components to 
be of comparable masses while in our calciilation fre- 
quencies are determined under assumption that the 
secondary is much less massive than the primary (cf. 
also Sundelius et al. 1997). As another example, a 
16-hr periodicity in the X-ray signal from the Seyfert 
galaxy IRAS 18325-5926 was described by Iwasawa et 
al. (1998). 

Let us suppose that the disk has only a weak dy- 
namical influence on the star and that stellar tides 
are negligible (the star is assumed much smaller than 
the typical curvature radius of the field around) as 
well as gravitational radiation and self-gravity both 
of the disk and of the star itself. Then the worldline 
of the star is very close to a geodesic in a pure gravi- 
tational field of the central rotating black hole. This 
is described by the Kerr metric which reads (Misner, 
Thorne & Wheeler 1973, p. 878), in Boyer-Lindquist 
spheroidal coordinates (t, r, 6, (p), in geometrized units 
(in which c = G = 1, c being the speed of light in vac- 
uum and G the gravitational constant) and with the 
( — h- 1— 1-) signature of the metric tensor {g^ii/), 

= _^dt2 + ^ sin2^(d(/)-a;Kdt)2 

+■|dr2+Sd^^ (1) 

where M and a denote mass and specific rotational 
angular momentum of the source and 

A = r^-2Mr-^a^ S = cos^ 6*, (2) 

A = {r^ + a^f - Aa^ sin^ 9, luk = 2Mar/A. (3) 

It has been demonstrated (Syer, Clarke, & Rees 
1991; Vokrouhlicky & Karas 1993, 1998) that un- 
der the above-described circumstances the stellar or- 
bit undergoes three secular changes: a decrease of 
the semi-major axis (the star spirals towards the cen- 
tre due to the loss of energy in collisions with the 



disk), circularization (the orbit becomes spherical, 
r = const), and a decrease of the amplitude of preces- 
sion of the orbital plane about the equatorial plane of 
the centre (the orbit gradually declines into the disk 
plane). Since the time scale for circularization is typ- 
ically found shorter than, or of the same order as, the 
time scale necessary to drag the orbit into the disk, 
one can expect that at late stages of evolution of the 
hole-disk-star system the star follows a nearly equa- 
torial spherical geodesic. This stage is also the one in 
which the star stays for a relatively long time. 

In the next section we discuss relevant properties 
of a precessing orbit in the Kerr spacetime. Then, in 
Sect. 3, we derive simple relations which appear to 
be appropriate for practical study of the relevant ob- 
jects. Equations further simplify if the object orbits 
not very close to the black hole, as discussed in Sect. 
4. Our method employs also the observed emission- 
line profiles which were vastly studied in recent litera- 
ture; we summarize the relevant formulae and results 
in Appendix A. 

2. The Orbiting Star as a "Pharaoh's Fan" 

For a general spherical geodesic in the Kerr space- 
time, the frequencies of the azimuthal and latitudi- 
nal motion are given by rather complicated formu- 
las involving elliptic integrals (Wilkins 1972; Karas 
& Vokrouhlicky 1994). However, the formulas sim- 
plify considerably in the case of a nearly equatorial 
geodesic. The azimuthal angular velocity lo = d(j)/dt 
(with respect to an observer at rest at radial infinity) 
can be approximated by that of the equatorial circular 
geodesic (Bardeen, Press, & Teukolsky 1972), 

LU± = {a + l/y±)-\ (4) 

where y± = y{u)±) = ±y/M/r^ are the correspond- 
ing values of the "reduced frequency" , in general de- 
fined by y = w/(l — aw) (de Felice 1994), and the 
upper /lower sign corresponds to prograde/retrograde 
orbit.® The angular frequency of (small) harmonic 
latitudinal oscillations about the equatorial plane is 
given, for a spherical orbit with general (but of course 
steady) radial component of acceleration, by 

n^={u'/ry{Au;' + 2yl[a-ir^+a^)u;]^}, (5) 



'Wc allow for both (±) cases for completeness, but solely the 
progradc trajectory can be considered: the accretion disk is 
more likely to be corotating with the central black hole and it 
has been shown that the interaction with the disk makes also 
the star corotate eventually. 
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where 



1 - 2MrS-^ (1 - aw sin^ ef 



(6) 



is the time component of the corresponding four- 
velocity, and Lo = const. The expression (5) was ob- 
tained by de FeUcc & Usscgho-Tomasset (1996), in 
studying the behaviour of a gcdanken apparatus bap- 
tized "Pharaoh's fan" (dc Felice & Usseglio-Tomasset 
1992), as a result of the analysis of the equation of 
relative deviation; for a simpler derivation by pertur- 
bation of the equatorial circular orbit, see Semerak & 
de Fehce (1997). 

We will briefly describe the "Pharaoh's fan" be- 
cause the equations describing its behaviour arc rele- 
vant also for the present work. The device consists 
of a monopole test particle in a non-friction short 
narrow pipe which constrains the particle's motion 
in two spatial dimensions but leaves it free in the re- 
maining direction. It was considered with the pur- 
pose to provide a space traveller with a measurement 
which, together with several other quasi-local experi- 
ments, yields a complete and non-ambiguous informa- 
tion about spacetime and about the orbit (Semerak 
& de Felice 1997). When the observer is moving on a 
(generally non-Keplerian) equatorial circular trajec- 
tory in the Kerr field and has the pipe turned to lat- 
itudinal direction, symmetrically to the orbital (i.e. 
equatorial) plane, the centre of the pipe is the par- 
ticle's stable equilibrium position. If released after 
being perturbed slightly off the equilibrium, the par- 
ticle harmonically oscillates ( "the Fan waves" ) about 
the equatorial plane, with the proper frequency 
This cannot be made vanish by any choice of the an- 
gular velocity w, since, in Kerr spacetime, it is not 
possible to orbit freely on circular orbits at a constant 
non-equatorial latitude (de Felice 1979). Except for 
the case a = (the Schwarzschild spacetime around a 
nonrotating centre) , 17^ is not equal to the square of 
the proper azimuthal frequency, (w*w)^, as a result of 
the "dragging" effect of mass currents existing within 
the rotating source (Lense & Thirring 1918; Wilkins 
1972). 

While the problem itself of the astronaut's ori- 
entability in black-hole fields is yet of a limited prac- 
tical usefulness, a question was raised (de Felice & 
Usseglio-Tomasset 1992; Semerak & de Felice 1997) 
whether some kind of correspondence could be estab- 



lished between the local experiments and the data 
which can be measured "at infinity". In the case of 
the AGNs or stellar-size X-ray sources, such a link 
is provided by the orbiting star as discussed above 
since it plays the role of the "Pharaoh's fan" ; in fact, 
the frequency of disk crossing would just be twice the 
Fan's frequency |f2|. (Since the star is free also in the 
radial direction, no pipe is necessary to constrain it.) 
We must also realize that what can in principle be 
inferred from (photometric) observations is not the 
proper frequency \Vt\, but rather the value measured 
by a distant observer, given by \^oo\ = i,Prom 
eqs. (4)-(5) the explicit expression is 



4(1 



4at/±-F3aV''^)- 



(7) 



3. The Set of Equations for Determining 
a and r 

A different kind of information can be obtained 
from the analysis made by Fanton, Calvani, de Fe- 
lice, & Cadez (1997, Sect. 7) relating the quantities 
which are observable in the integrated spectrum of 
an accretion disk to the parameters of the hole-disk 
system. A stationary disk produces the well-known 
double-horn line profile which is presumably modu- 
lated by the star-disk interaction at a certain radius 
(see Appendix for details). Approximating the star 
as an emitting point-like source on a circular equato- 
rial geodesic, the observed frequency shift g of each 
emitted photon can be written in terms of its direc- 
tion cosinus (azimuthal component of the unit vec- 
tor along the direction of emission of a given photon, 
measured in the emitter's local rest frame) as 



l_2M/r + t/± (a + VAe^)] , (8) 



where 



1 - Md/r + 2a?y± (9) 

from (5) and (4). One of the most important at- 
tributes of a spectral line is its width; this arises from 
the different frequency shifts g carried by the photons 
which reach the observer at infinity. Since g depends 
on the emission cosinus e^, the maximum broadening 
of the line which could be measured as a result of an 
integration over one entire orbit, is a measure of the 
variation 5g corresponding to the maximum range of 
variability of compatible with detection at infinity. 
If we call the latter i, we have from (8) 



{5g) = (^e^)u*a;±\/A 



(10) 
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and thus 



2 2 

= w I r 



1 - 2M/r + a^/r'^ 
1 - 3M/r + 2ay± 



(11) 



It is clear that Se^ can be at most 2 but, in rcaUstic 
situations, it varies significantly with the inclination 
angle of the black hole-disk system with respect 
to the line of sight. (It also depends, though only 
weakly, on the rotational parameter a and the radius 
of emission r.) From a numerical ray-tracing analysis 
(Fanton 1997, private communication) it is found that 
Se^ ranges from ((5e^)niin ^ 0.4 to ((5e^)inax ^ 2 as 
6o goes from ~ 0° to ~ 90°, hence one can fix as 
observables the extreme values of S, (^max = 2.5 Sg and 
Smin = 0.5 Sg, (corresponding to a line of sight nearly 
polar in the former case and nearly equatorial in the 
latter one). If the line width 5g is measured, then 
formulas (4), (7) and (11) provide a closed system of 
ordinary equations which yields the parameters M, a 
and r in terms of the observable quantities u}±, |0oo| 
and 6. 

These equations can be solved for a, and y±, for 
example: 



a 

„2 



= w. 



± y± ' 



5loI' 



Auj±y± + 

where y± is determined by quartic equation 



Myl - Byl + Cyl-Dy±+E = 



with 



B = 

C = 

D = 

E = 



— (5'a;^ + 23a;i-20^), 
215^0;^ + 76a;| + 7(5^0^ + 30^ 

ui± ^ =^ =^ =^ °° 

+11^1^1, 
5a;^(452 + 3) 

4 



(^2 + 1) n 



^ 1 

ooj ' 



+(a;|-n2^)(^2j^2^-7a;^)-f2t, 



(12) 
(13) 

(14) 

(15) 
(16) 

(17) 
(18) 



Evidently it would be rather cumbersome to solve eq. 
(14) analytically and one better finds the result nu- 
merically for each given set of data w±, \^oo\ and \5\. 

Let us note that in the Schwarzschild case, a = 0, 
eqs. (4) and (7) reduce to 



a/M 
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40.4 


0.85 


5 


0.0831 


0.0735 


0.925 


5.3 




10 


0.0308 


0.0294 


0.868 


10.3 




20 


0.0110 


0.0108 


0.855 


20.4 




40 


0.0039 


0.0039 


0.851 


40.4 


1.00 


5 


0.0821 


0.0716 


1.200 


5.3 




10 


0.0306 


0.0291 


1.033 


10.3 




20 


0.0110 


0.0108 


1.008 


20.3 




40 


0.0039 


0.0039 


1.002 


40.4 



Table 1: Illustration of accuracy of approximative re- 
lations; a/M and r/M arc estimates based on cqs. 
(25)-(26). Notice that acceptable precision is reached 
for r ^ lOM. 



and eq. (11) reads 



1 - 2r^ujl 
1 - Sr'^wl 



(20) 



The physical solution of this equation is 



r2 = (4w^)-i 36^ + 1- V(3(52 + 1)2 - 86^ , (21) 
M following then from (19) as M = r^oj^.. 

4. A Simple Explicit Solution for r^';$>a^ 

Due to processes of tidal disruption close to the 
horizon (Frank & Rees 1976; Marck, Lioure, & Bonaz- 
zola 1996; Diener et al. 1997) and in general due to 
rather violent conditions near the inner edge of the 
disk'' where most of the radiation is produced, the 
star probably cannot orbit, in the manner required 
above, on radii of the order of the centre's gravita- 
tional radius (~ M). On the other hand, the radius 



(19) 



^This is standardly considered somewhere near the marginally 
stable and marginally bound circular equatorial geodesic orbits 
which may lie, in terms of the r-coordinate, very close to the 
horizon (r ~ M) when the black hole rotates rapidly (a ~ M). 
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must not bc' too large' since this would prevent us to 
recognize any relativistic effect and even to cover a 
sufficient number of periods in the observations. It is 
then relevant to suppose that the star orbits at some 
r ^ lOM, and typically r ~ 30M. 

One can verify that some terms in eqs. (4), (7) 
and (11) are then very small when compared with 
the rest. Ignoring these small terms, one arrives at a 
simplified set of relations which keeps an acceptable 
precision of the resulting parameters without need to 
solve the fourth-order equation (14). Exact equations 
can only then be used to improve the estimates. It 
is suitable to choose some particular approximation 
according to what precision one desires at each of the 
quantities M, a, r, and also according to precision of 
each of the input data \oJ±\, \^oo\, |^|- For instance, 
since ay± < 1/30, we can take 



w± = y±{l - ay±), y± = oj±{l + au±) (22) 

(for the determination of M from the known r, or vice 
versa, it is even sufficient to neglect ay± altogether 

2/2 

and use just Schwarzschild equation). Also 2M/r ^ 
but the other relevant bounds are rather weak, 
namely 

a^/r^ _ ay± _ a 1 
/Mr 5 ' 



(23) 



ay± M/ r 

so one must be careful when making neglections in 
eqs. (7) and (11); in particular, neglecting the term 
a^/r^ in eq. (7), a very simple formula for a is ob- 
tained just by combination with eq. (4), 



a 



1 



(24) 



but this turns out to yield inaccurate estimates. 

Let us present an example of a simple and satis- 
factory approximation. (We denote the approxima- 
tive values by a tilde.) Suppose the observed signal 
is produced by a star on a prograde orbit {plus sign 
in relevant relations). Setting the fraction in eq. (11) 
equal to unity, one obtains an approximative value for 
the radius, 



= S/uj+; 



(25) 



within the considered range of a and r, the error is 
6.5% at maximum (for a = and r = lOM) but 
typically only 1.3% (corresponding to a = 0.8M, r = 



30M ) . Substituting f and cj+ in place of r and y+ in 
eq. (7), one approximates a by 



3u+ 



25 



^4<52-3(l-J12^K) 



M is then given by eq. (4): 



5^ 



+ a;+(l + aw+)^ 



(26) 



(27) 



The most desired information, namely the values of 
a/M and f/M, follow by obvious combinations from 
the deduced a, r and M. The results of this approx- 
imation are illustrated in Table 1 for several typical 
a and r, together with the respective values of the 
observable frequencies. 

Before using the above formulas for some partic- 
ular set of real data, one must realize that they are 
written in gcomctrizcd imits. We remark that the rel- 
evant quantities in physical units are obtained by the 
following conversions: 



M 



M, 



^phys 



phys 



ca. 



1.477 X lO^cm ' 



r-phys 



also 



7 phys 



„phys 



(28) 
(29) 



(30) 



To obtain frequency in physical units [Hz] (either from 
oj or Q), one uses the relation 



^phys _ 



27r 27r 



(4.771 X 10^cm-s-i)w. (31) 



In usual notation of relativistic astrophysics the (ge- 
ometrized) frequencies are scaled by M~^, as in 
Tab. 1. Numerical values thus obtained must be mul- 
tiplied by a factor 

— ^ = (3.231 X lO^Hz) (32) 
to get the frequency in [Hz] . 

5. Discussion 

In the present paper, we outlined general features 
of our method; an application to particular sources 
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will be discussed separately (work in progress). We 
will conclude by brief comments on several important 
points. 

First, although the method described above in 
principle yields all the relevant quantities, it might 
be combined advantageously with independent deter- 
mination of some of the parameters. In particular, 
knowing M, one can deduce a/M and r/M from eqs. 
(26) and (25). 

Observed emission-line profiles from relativistic ac- 
cretion disks around black holes have recently at- 
tracted great interest and they play a significant role 
also in our method. Apart from the well-known 
double- horn feature (Laor 1991), characteristics of the 
lines depend on rather uncertain properties of accre- 
tion flows, e.g. on advection velocity (Fukue & Ohna 
1997), limb-darkening law (Rybicki & Bromley 1997) 
and shape of the disk (Pariev & Bromley 1997). One 
therefore needs to investigate a broad range of mod- 
els to determine widths, required in our model, and to 
reject unacceptable profiles. We present some results 
in this direction in the Appendix. 

Crucial point in our considerations is the (quasi)- 
periodic variability of the black-hole-disk source. Pe- 
riodic modulation of the accretion flow is the most 
likely signature for a stellar companion close to the 
central black hole (Podsiadlowski & Rees 1994). Each 
star's passage through the accretion disk pulls some 
amount of gaseous material out of the disk (Zurek, 
Siemiginowska, & Colgate 1994). This material tem- 
porarily covers the innermost region of the disk and 
modulates the observed radiation (both contimmm 
and line).^ It can thus be anticipated that a variable 
signal is recognized in the X-ray band, however, de- 
tails of the mechanism remain rather imcertain. Note 
that also the star may be affected by the collisions 
and tidal forces considerably — it may lose its outer 
layers or even become disrupted (Frank & Rees 1976; 
Marck, Lioure, & Bonazzola 1996; Diener et al. 1997; 
Loeb & Ulmer 1997). A stripped stellar core is an 
intense source of (ir) radiation which may survive for 
many orbital periods (Rees 1998). In any case, we 
assumed here that the star is only weakly affected by 
the disk, i.e. that its orbit may well be approximated 
by a nearly equatorial spherical geodesic. 

Radiation from the presumed accretion disks can 



Syer & Clarke (1995) discussed a response of the disk to a 
body moving completely inside and determined conditions un- 
der which a gap can develop and survive (see also glechta 1998). 



be modulated by their precession, vibrations and 
various instabilities. It is also very well possible 
that the disks are covered by numerous irregulari- 
ties (phenomenologically designated as bright spots), 
contributing to the featureless short-term X-ray vari- 
ability which is indeed observed (Abramowicz et al. 
1991; Mangalam & Wiita 1993). Then our scheme 
must anticipate that well-separated, profound irreg- 
ularities develop on the disk surface (due to stellar 
passage) and survive several orbital periods. 

Although no AGN or black-hole X-ray binary with 
clear periodic variability is currently known, the cases 
of rapidly variable Ka line (Iwasawa et al. 1996; 
Yaqoob et al. 1996), and that of quasi-periodic AGNs 
(Papadakis & Lawrence 1993;^ Lehto & Valtonen 
1996; Stothers & Sillanpaa 1997; Valtonen & Lehto 
1997) and X-ray sources (Callananet al. 1992; Pavlen- 
ko et al. 1996; Belloni et al. 1997; Steiman- Cameron 
& Scargle 1997; Iwasawa et al. 1998) suggest that 
some periodic components are present in the signal. 
Thus we conclude with M. J. Rees (1998): "There is a 
reaJ chance that someday observers will End evidence 
that an AGN is being modulated by an orbiting star, 
which would act as a test particle whose orbital pre- 
cession would probe the metric in the domain where 
the distinctive features of the Kerr geometry should 
show up clearly." 
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A. Appendix 

This Appendix summarizes effects which influence 
formation of observed line profiles. We have system- 
atically studied a large set of synthetic emission-line 
profiles from an accretion disk around a rotating black 
hole which are needed for comparisons with observa- 
tional data in the above-described method. Various 
properties of observed radiation from a source near 
a black hole have been discussed by many authors 
using different approximations, starting with classi- 
cal papers by Cunningham (1976), Laor (1991), and 
more recently by, e.g., Fanton et al. (1997). Our aim 
was twofold: (i) to determine the range of redshift fac- 
tor affecting the predicted line profiles under different 
situations; (ii) to see how sensitive the results are to 
uncertainties in current models of the disk emission. 
We considered regions with different radial distance 
from the central black hole and took into account 
also non-negligible radial motion of gaseous material, 
different local profiles of emission lines and the disk 
shape. Typical results are illustrated in Figures Al- 
A4. 

The model of the; disk is characterized by radiation 
intensity Jem emitted locally from its surface in the 
equatorial plane at frequency Vem, 

) =F(rer„)¥'l(l'em)¥'2(t?), (Al) 

where -F(rcm) is total radiation flux at the disk sur- 
face, (pi{i'em) is the emissivity profile in frequency, 
ip2{'&) is the limb- darkening law; i? denotes the angle 
between the ray and direction normal to the disk in 
LDF. We have classified different models according to 
following parameters: 

1. Local radiation flux, f (rem)- We examined: 
(i) power- law dependence i^(rem) oc r^m^^^ (/3 = 
const); (ii) standard Novikov & Thorne (1973) 
model. 

2. Emissivity profile in frequency, (^i(fc'em)- We ex- 
amined two cases: (i) Gaussian profiles </7i(fem) 
oc exp[— e(z/om ^ 1)]^ (e ~ const). The frequency 
profile is normalized to maximum at unit local 
frequency [thus the term (fem — 1)]. (ii) Asym- 
metric profiles which account for self-irradiation 
of the disk (this topic was discussed by several 
authors but effects of irradiation of the disk are 
still an open problem; cf. Petrucci & Henri 1997; 
Dabrowski et al. 1997). 
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3. Angular dependence of the local emissivity (limb- 
darkening law), <p2(^)- We examined: (i) 1^2(1^) cx: 
1 + e/i; jj, = cos 9, e = const (e = corresponds 
to locally isotropic emission). (ii) (/?2(^) c< 
Atlog(l+l/At) (Ghisellini, Haardt, & Matt 1994). 

Ray-tracing of the light trajectories was performed 
in the Kerr geometry with corresponding redshift 
function g being given by 



PaV 

-g'^^^.m + 9""flrPr/Pt , 



and 



cos ^ = — 



Pan" 

PaV" 



(A2) 



(A3) 



Here, and 77" denote respectively the four-momenta 
of the photon and of the emitting material in the disk, 
and analogously p" and fj" that of the photon and of 
the observer at r ^ 00; n" denotes a unit space-like 
vector perpendicular to the disk surface. 

The resulting synthetic lines are characterized by 

observed radiation flux (count rate) as a function of 
energy. Predicted profiles are shown in figures A1-A4 
(background is subtracted and flux is normalized to 
maximum) . Characteristics of the line proflles can be 
represented in terms of the following parameters: 

1. Normalized centroid energy E"^ (Matt, Perola & 
Piro 1991) which is determined by energy dis- 
tribution N{E) of the line photons: 

E'' = N-^ j N{E)EdE, N = j N{E)dE; 

(A4) 

E"^ characterizes the shift of the line energy (no- 
tice that E"^ is not the standard energy centroid 
which is defined through the continuum flux). 

2. Ratio of line widths (of the observed width to 
the width in LDF) 6a. Here, geometrical line 
width <t{E'') is defined by 

a'^iE'') = AT-i J N{E){E - E^fdE . (A5) 

Broadening of the line corresponds to Sa > 1. 

3. Normalized (to the unit continuum flux) equiv- 
alent width 



W = j N{E)EdE. 



4. Energy F"^ of the line center, given by 



J 

Jo 



J 

JF' 



N{E)EAE= / N{E)EAE. (A7) 



5. Energy of the maximum count rate, F™. 

In Fig. Al, the emitted line profile (in the local 
frame corotating with the disk material) is indicated 
by shading: here, it is the Gaussian line originating in 
a thin Keplerian disk. Photon energy is normalized 
to unity and indicated on the horizontal axis. Locally 
emitted radiation flux is isotropic (e = 0) and de- 
creases as oc r~^/^. The thick solid curve corresponds 
to the observed profile from the whole disk which ex- 
tends between rjn = QAd < r < 25M = rout- Contri- 
butions from three equally wide regions between rin 
and Tout are shown by thin solid curve (inner part of 
the disk), dot-dashed curve (central part), and the 
dashed curve (outer part). In Fig. A2, parameters 
are the same except for local velocity of the mate- 
rial which has advective component Vr = = 0.4w+; 
v+ is the local corotating Keplerian orbital velocity, 
Vr and are three-velocity components in the frame 
corotating with the disk material. (Numerical value 
of v,j,/vr is model-dependent; here, it has been cho- 
sen for the sake of comparison with Fukue & Ohna 
1997; Watanabc & Fukue 1996 suggested an alterna- 
tive scheme of an advective corona.) Profiles given in 
Fig. A3 are constructed as in Fig. Al but for locally 
anisotropic radiation (e = 2), and similarly Fig. A4 
illustrates profiles from the disk with both radial fiow 
and locally anisotropic radiation. 

Parameters characterizing geometrical shapes end 
energy shift of the profiles from all four figures are 
summarized in Tab. 2. All the line profiles in this 
example have F§ — Eq ~ 1, VFo = 105 in the lo- 
cal disk frame (subscript "0"). Much more extensive 
set of predicted profiles and corresponding tables are 
available in the electronic form from the authors. 



(A6) 
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a/M=0 a/M= 0.5 a/M=^ 




Energy 

Fig. Al. — Synthetic emission-line profiles from a thin Keplerian disk with locally isotropic radiation and for 
three values of the black-hole dimensionless angular momentum: a/M — (left column), a/AI = 0.5 (center), and 
a/M — 1 (right). Observer's inclination 9o acquires values between 9o = 0° (pole-on view) and 90° (edge-on). See 
the text for a detailed description and Tab. 2 for characteristics of the profiles. 
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Fig. A2. — As in Fig. Al but for a disk with radial component of the flow velocity comparable to the azimuthal 
component. 
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Fig. A3. — As in Fig. Al but for locally anisotropic radiation. 
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Fig. A4. — As in Fig. A2 but for locally anisotropic radiation. 
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TABLE 2 

Characteristics of the line profiles. The first column refers to the corresponding figure 



Total 



Inner 



Central 



Outer 



617,7 
77e.^ 
1329.6 
2446.1 
2728.4 
596. C 
771.7 
1325.9 
2429.5 
2654.5 
577.7 
769.6 
1322.5 
2413.5 
261Q.6 
576.4 
543.2 
7QQ.7 
1375.3 
1873. S 
568.3 
544.2 
713.1 
1414.2 
1887. C 
561.8 
546.5 
727.9 
1457. C 
1880. £ 
1041. C 
1306.5 
2115.7 
3428.8 
2933.7 
1007.2 
130Q.6 
2110.6 
3403.1 
2838.7 
978.5 
1298.4 
2105.8 
3379.3 
2790.9 
983.8 
912.8 
1109. C 
1918.4 
2090.6 
970.8 
914.6 
1128.4 
1972. £ 
2093.9 
960.2 
918.6 
1151.5 
2032.3 
2075.8 



W 

470. 

941. 

953. 
1433. 
2191. 

466. 

934. 

949. 

2168. 

463. 

925. 

945. 
1395. 
2133. 

386. 

736. 

852. 

851. 

822. 

385. 

731. 

839. 

829. 

845. 

385. 

726. 

823. 

806. 

884. 

465. 

905. 

915. 
1430. 
2508. 

462. 

898. 

911. 
1407. 
2489. 

459. 

890. 

908. 
1390. 
2440. 

384. 

718. 

826. 

846. 

893. 

383. 

714. 

814. 

825. 

934. 

383. 

710. 

800. 

802. 

992. 



32G.6 
525. C 
1Q75.6 
230G.6 
2646.2 
313. -3 

1084.: 
2298. C 
..19 2597.7 
302.3 
541.7 
1090.8 
228,1 
2G05 
.94 310 



348.3 
525.8 
1157.3 
1717.9 
299. C 
358.8 
557.6 
1227.8 
1697.7 



1.04 



545.5 
883. -3 
171G.3 
32G8.: 
..19 2872. -3 
525.8 
898.: 
1730.6 
3253. C 
..19 2807.8 
508.: 
913.5 
1741.^ 
..IG 3227. C 
2819. -3 



W 

530. 
1094. 
1008. 
13G2. 
2334. 

517. 
1062. 

984, 
1309. 
2187, 

502. 
1029. 

9G5. 
1 27,T 
2014, 



920, 
791. 
41G. 
83G. 
931. 
873. 
831. 

820. 

885. 

833. 

887. 

525. 
1041. 

9G4. 
1359. 
2740. 

513. 
1013. 

94,1 
1305, 
2549. 

499. 

985. 



.94 



527.5 
572..- 
79r,,1 

15.17..'^ 

1908.7 
519.1 
S8G.: 
837.: 

1628.7 

I89r, 

509 
G04.3 
88G.9 
1725. -3 
1837.: 



813, 
899. 
872. 

412 
79G. 
857. 
834. 
1057. 



L.12 1.03 



233.8 
530.8 
904.: 
1815.^ 

2033. * 
232.: 
530.^ 
899. -3 

1809.9 

2034. : 
231 
529.^ 
898.5 

1 802 f 
2029, t 
216,8 
309.8 
460.1 
915.^ 
997.8 

217. r 

312.9 
467.1 
917.5 
1003.^ 
217.1 
316.1 



398.^ 
901.8 
1443.1 

2530.: 
2226.8 

396. e 
901.1 
14,15 : 

2524,8 
2228.^ 
394,5 
899.^ 
1433. e 
2516.9 
2224.: 



711 8 
1294.1 
1102.6 
372.8 
530. C 
744.1 
1294.; 
1106.'; 
373. C 
535.1 
756.: 
1299.: 



353. 

772. 

900. 
1001. 
1043. 

350. 

768. 

89G. 

989. 
1023. 

349. 

766. 

892. 

983. 
1012. 

271. 

591. 

708. 

685. 

690. 

270. 

592. 

703. 

693. 

692. 

2 70. 

593. 



697. 

35.1. 

742. 

860. 

96G. 
1028. 

350. 

738. 

8 57. 

953. 
1007. 

349. 

73G. 

854. 

948. 

99G. 

271. 



6 54. 
6 73. 
270. 
577. 
677. 
664. 
677. 
2C9. 
579. 
67G. 
672. 
684. 



209.9 
4G3.4 
757.9 
1509.7 
1661.7 
209.7 
462.4 
757.1 
1502.1 
1G54.4 
209.4 
461. C 
757.3 
1 497.2 
1G4E.2 
198.6 
273.7 
384.7 
G7G.8 
G98.C 
198.7 
274.5 
38G.1 
G81,8 
702.8 
198.6 
275.2 



359 

788 
1207. 
2102 
1819. 

359 

78G 
1 20,':; 

2092 

1811 
358 
784 

1206 

2085 

1801 
342.1 
4G3.5 
G12.4 
948.5 
7G4.7 
342.3 
4G4.7 
G14.C 
954.4 
770.3 
342.2 
4G5.8 
G17.5 
9G0.3 
777.3 



307. 
731. 
830. 
833. 
865. 
306. 
731. 
827. 
832. 
8G3. 
306. 
728. 
823. 
813 
863. 
238. 
546. 
668. 
724. 
794. 
238. 
550. 
G73. 
722. 
793. 
238. 
551. 
GO G74. 
8G 722. 
54 788. 
69 307. 
45 705. 
OG 793. 
8G 794. 
39 850. 
48 306. 
75 705. 
90 791. 
4G 794. 
31 848. 
8G 306. 
13 702. 
06 787. 
92 794. 



238. 
535. 
G44. 
691. 
781. 
238. 
539. 
G49. 
691. 
780. 
238. 
540. 
G51. 
G92. 
776. 



NOTE. — a* . . . dimensionless black-hole angular momentumi Aq. . . observer inclination; Aa. . . Ratio of the widths o£ the synthetic line (observed versus LDF); W. . . Normalized equivalent width of the synthetic line; 
. . Energy q£ the line center; . . Centroid energy; F™. . . Energy q£ the maximu m count rate (all quantities are nonualizedi aubacript "0" correaponda to the LDF). 



